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Variational Monte Carlo and Green's function Monte Carlo are powerful tools for cal- 
culations of properties of light nuclei using realistic two-nucleon (NN) and three-nucleon 
(NNN) potentials. Recently the GFMC method has been extended to multiple states 
{Sj | with the same quantum numbers. The combination of the Argonne Vi$ two-nucleon and 
Illinois-2 three-nucleon potentials gives a good prediction of many energies of nuclei up 
to 12 C. A number of other recent results are presented: comparison of binding energies 
with those obtained by the no-core shell model; the incompatibility of modern nuclear 
Hamiltonians with a bound tetra-neutron; difficulties in computing RMS radii of very 
weakly bound nuclei, such as 6 He; center-of-mass effects on spectroscopic factors; and the 
possible use of an artificial external well in calculations of neutron-rich isotopes. 
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1. INTRODUCTION 



In the last decade, the Green's function Monte Carlo (GFMC) method has been de- 
veloped into a powerful tool for calculations of light nuclei (so far up to A = 12) using 
realistic two-nucleon (NN) and three-nucleon (NNN) potentials. GFMC starts with a 
trial wave function that is obtained via variational Monte Carlo (VMC) and projects out 
excited state contamination to, in principle, obtain the true ground-state wave function 
for the given Hamiltonian. In practice the method obtains ground and low-lying excited 
state energies with an accuracy of 1-2%. A review of the nuclear VMC and GFMC meth- 
ods up to A = 8 may be found in Ref. [H]; A = 9,10 results are in Ref. [|2j. By using 
the Argonne vi$ NN potential (AV18) and including two- and three-pion exchange NNN 
potentials, a series of model Hamiltonians (the Illinois models) were constructed that 
give a good reproduction of energies for A = 3 to 12. Recently GFMC has been extended 
to the calculation of multiple excited states with the same quantum numbers [Ej. 

This paper reviews some of these results and compares GFMC energies to No-core Shell 
Model (NCSM) results [|3] for several cases. A recent study showing that an experimental 
claim of a bound tetraneutron is very unlikely to be valid [ E] is also reprised. Finally 
some on-going work involving spectroscopic factors, neutron-rich oxygen isotopes, and 
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difficulties in GFMC computation of RMS radii is presented. 

2. HAMILTONIANS 

Our Hamiltonian includes a nonrelativistic one-body kinetic energy, the Argonne Dig 
(AV18) two-nucleon potential [Ej and various three-nucleon potentials, 

h = U-^) + E^+ E (i) 

i LUl i<j i<j<k 

The difference between proton and neutron masses is included in our calculations, but 
ignored above. The Argonne v-]_& potential is one of a number of accurate AW potential 
models developed since 1990. It can be written as a sum of electromagnetic and one- 
pion-exchange terms and a shorter-range phenomenological part. The electromagnetic 
terms include one- and two-photon-exchange Coulomb interactions, vacuum polarization, 
Darwin-Foldy, and magnetic moment terms, with appropriate proton and neutron form 
factors. 

The one-pion-exchange part contains the normal Yukawa and tensor functions with a 
short-range cutoff. This and the remaining phenomenological part of the potential can 
be written as a sum of eighteen operators, which is where the name V\g comes from. The 
first fourteen are charge-independent, and include spin-spin, tensor, L- S, and quadratic- L 
terms, each with a dependence on isospin. The last four operators break charge indepen- 
dence. The radial forms associated with each operator are determined by fitting AW 
scattering data. The potential was fit directly to the Nijmegen AW scattering data base [ 
|H], which contains 1787 pp and 2514 np data in the range 0—350 MeV, with a \ 2 P er datum 
of 1.09. It was also fit to the nn scattering length measured in c?(7r~,7)nn experiments 
and the deuteron binding energy. 
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Figure 1. Terms in the Illinois three-nucleon potentials. 



The Urbana series of three-nucleon potentials were developed to compute properties 
of A = 3,4 nuclei and nuclear matter; the current (ninth) model is designated UIX 
These potentials are written as sums of two-pion-exchange with intermediate excitation 
of an isobar (left panel of Fig. and shorter-range phenomenological terms. The two- 
pion-exchange term is that of the original Fujita-Miyazawa model and contains both 
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spin (tensor) and isospin dependence. The shorter-range phenomenological term is purely 
central and repulsive. Our recent work has shown the need for additional binding for 
p-shell nuclei and for further increased binding as \N — Z\ increases. This led to the 
development of the Illinois models [Ej which, in addition to the Urbana terms, contain 
the two-pion s-wave scattering term (second panel of Fig. and three-pion exchange 
ring terms (last two panels of Fig. [TJ. The latter can involve the excitation of one or two 
sequential isobars, so that each energy denominator contains only one A mass. 

In light nuclei we find that the three-nucleon potential contributes only 2-9% (increasing 
with A) of the total potential energy. However, due to the large cancellation of potential 
and kinetic energy, this amounts to 15-50% (increasing mostly with N — Z) of the binding 
energy. We expect a similar ratio for the four-body potential, which implies that it 
contributes only a few percent of the binding energy; such contributions are close to our 
computational accuracy and would be absorbed in fit of the NNN potential parameters. 

3. QUANTUM MONTE CARLO METHODS 

3.1. Variational Monte Carlo 

Variational Monte Carlo finds an upper bound, Et, to an eigenenergy of H by evaluating 
the expectation value of if in a trial wave function, iff?. The parameters in are varied 
to minimize E?, and the lowest value is taken as the approximate energy. Over the years, 
we have developed rather sophisticated for light nuclei [|TTJ[T2|. A good variational 
trial function has the form 



|*T> 



i+ E u™ 1 



i<j<k 



i<j 



■ (2) 



The Uij and U™ 1 are non-commuting two- and three-nucleon correlation operators (the 
most important being the tensor-isospin correlation corresponding to the pion-exchange 
potential); S indicates a symmetric sum over all possible orderings; and tyj is a fully 
antisymmetric Jastrow wave function which determines the quantum numbers of the state 
being computed. 

The Jastrow wave function, j, for p-shell nuclei starts with a sum over independent- 
particle terms, that have 4 nucleons in an a-like core and (A — 4) nucleons in p-shell 
orbitals. These orbitals are coupled in a LS[n] basis to obtain the desired JM value 
of a given state, where n specifies the spatial symmetry [n] of the angular momentum 
coupling of the p-shell nucleons. Different possible L5[ra] combinations lead to multiple 
components in the Jastrow wave function. This independent-particle basis is acted on by 
products of central pair and triplet correlation functions: 

m = m[yi #j[ n unj)} e (fowl n fs P s[n \r kl ) 

\ i<j<k i<j<4 LS[n] k<i<l<A 



X 



II tf/ W (^)ll*A(^[n]JMTT 3 ) 123 4 : 5...A)) • (3) 



pp 

4<l<m<A 



The operator A indicates an antisymmetric sum over all possible partitions of the A 
particles into 4 s-shell and (A — 4) p-shell ones. The pair correlation for both particles 
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within the s-shell, / ss , is similar to that in an a-particle. The pair correlations for both 
particles in the p-shell, fp P s ^ n \ and for mixed pairs, fg P s ^ n \ are similar to f ss at short 
distance, but their long-range structure is adjusted to give appropriate clustering behavior, 
and they may vary with LS[n]. 

The single-particle wave functions $a are given by: 



<5> A (LS[n]JMTT 3 ) l2U:5 ... A ) = |$ 4 (00 00) 1234 [ II </>f W 0M] (4) 

4<1<A 

[ n . Y ^ai)] LMiM [ n . x4m s )] [ n ^)] TT ) ■ 



4</<A 1 ' 4</<A ) JM 4<1 < A 

The (f)p S ^ n \R a i) are p-wave solutions of a particle in an effective a — N potential that has 
Woods-Saxon and Coulomb parts. They are functions of the distance between the center 
of mass of the a core and nucleon I, and may vary with LS[n}. The depth, width, and 
surface thickness of the single-particle potential are additional variational parameters of 
the trial function. The overall wave function is translationally invariant, so there is no 
spurious center of mass motion. 

The @LS[n] mixing parameters of Eq. (J3J) are determined by a diagonalization procedure, 
in which matrix elements 

E TM = (^ T (fti)\H\9 T (fij)) , (5) 
N T4j = <* t (A)|*t(&)) , (6) 

are evaluated using trial functions \&t(A) = ^t(A = 1,/%^ = 0). Although the 
$ a{LS[ti\J MTT3) are orthogonal due to spatial symmetry, the pair and triplet corre- 
lations in make the different LS[n] components nonorthogonal, so a generalized eigen- 
value routine is necessary to carry out the diagonalization. 

3.2. Green's Function Monte Carlo 

GFMC projects out the lowest-energy ground state from the VMC \I/t by using 

tf(r) = e - {H - Eo)T m T ; (7) 
^0 = Km *(r). (8) 

T— >0O 

If sufficiently large r is reached, the eigenvalue Eq is calculated exactly while other ex- 
pectation values are generally calculated neglecting terms of order — \&t| 2 and higher. 
In contrast, the error in the variational energy, Et, is of order \^q — ^t| 2 , and other 
expectation values calculated with ty? have errors of order \ty — ^>t\- Here I present 
a simplified overview of nuclear GFMC; a rather complete discussion may be found in [ 

HUE!. 

Introducing a small time step, At, r = nAr, gives (typically Ar = 0.5 GeV -1 ) 



e -(H-E )Ar 



^ T = G n ^ T . (9) 



where G is the short-time Green's function. The \P(t) is represented by a vector function 
of R, and the Green's function, G aj g(R', R) is a matrix function of R' and R in spin-isospin 
space (labeled by the subscripts a,j3), defined as 

G aP (R', R) = (R, a\e- (H - Eo)AT \R, /?). (10) 
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Omitting spin-isospin indices for brevity, \E'(R n ,r) is given by 

*(R n ,r) = J G(R n ,R n _ 1 )---G(R 1 ,R )^ r (R ) C /R n - 1 ---rfRirfRo, (11) 

and the integration is done by Monte Carlo. We approximate the short-time propagator 
as a symmetrized product of exact two-body propagators and include the to first- 
order. In recent benchmark calculations [Ej of 4 He using an eight-operator NN potential, 
the GFMC energy had a statistical error of only 20 keV and agreed with the other best 
results to this accuracy (< 0.1%). 

For more than four nucleons, GFMC calculations suffer significantly from the well- 
known Fermion sign problem. This results in exponential growth of the statistical errors 
as one propagates to larger r, or as A is increased. For A > 8 the resulting limit on r 
is too small to allow convergence of the energy. In the last few years we have developed 
and extensively tested a constrained-path algorithm for nuclear GFMC [ H2] to solve this 
problem. In this method configurations with small or negative ^(ry-^T are discarded such 
that the average over all discarded configurations of ty(ry -^t is 0. This means that, if ^/t 
were the true eigenstate, the discarded configurations would contribute nothing but noise 
to (H). This constrained propagation completely controls the growth of the statistical 
errors and in most cases produces a result that is statistically the same as unconstrained 
propagation (the accuracy of the comparison may be limited by the statistical errors in the 
unconstrained result). However we have demonstrated some cases for which constrained 
propagation leads to a wrong result, and in fact for which the approximate (H) is not 
even an upper bound to the correct eigenvalue. In all cases the correct result can be 
obtained by making a few (10 to 20) unconstrained steps before evaluating the energy. 
Our calculations for A > 5 are now all made using constrained-path propagation with 10 
to 20 unconstrained steps. 

The number of spin-isospin components in ty? grows rapidly with the number of nu- 
cleons. Thus a calculation of a state in 8 Be involves about 30 times more floating-point 
operations than one for 6 Li, and 10 Be requires 50 times more than 8 Be. Calculations 
of the sort being described here are currently feasible up to only A = 12; for v4=10, 
these require ~8,000 processor hours on the NERSC IBM SP (Seaborg) running at 390 
MFLOPS/processor (10 16 operations) and ~150,000 processor hours for 12 C on the Los 
Alamos qsc computer running at 360 MFLOPS/processor (2 x 10 17 operations). 

Figure 121 compares the VMC and GFMC energies of several nuclei for the AV18+IL2 
Hamiltonian. We see that the variational wave functions for the s-shell nuclei are quite 
accurate-the GFMC improves our best VMC energy of 4 He by only 1.5 MeV or 5%. 
(These do not contain explicit correlations for the new terms in the Illinois NNN 
potential; the corresponding error for the AV18+UIX Hamiltonian is only 2%.) However, 
the p-shell variational wave functions are much less accurate; compared to the GFMC 
energies, they result in underbindings of 4.3 MeV (13%) for 6 Li to 25 MeV (38%) for 10 B. 
In fact the fail to give particle-stable energies for any of the A > 6 nuclei and give 
maximum binding energies for 8 Be. On the other hand, the excitation spectra from VMC 
and GFMC calculations are generally quite similar, except when there is a change in the 
dominant symmetry. 
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Figure 2. Comparison of ty? an d GFMC energies for several nuclei 



3.3. GFMC Evaluation of Excited States 

It is possible to treat at least a few excited states with the same quantum numbers 
using VMC and GFMC methods [H]. The VMC calculations have been described above, 
and essentially involve solving a generalized eigenvalue problem, Eqs. (0) and ©. The 
same basic method can be applied in GFMC calculations, though the implementation is 
slightly more involved. In this section, ^t,% represents the trial wave function for the i th 
state of specified (J w , T) and ^(r) is the GFMC wave function propagated from it. By 
construction {^T,i\^T,j) — for % ^ j. We would like to calculate the Hamiltonian and 
normalization matrix elements as a function of r: 



(^(r/2)| H |ttj(T/2)) 

iA) ~ 1^/2)11^/2)1 ' {l2) 
_ (^(r/2)|^-(r/2)) 

! ^ [T) 1*^/2)11^(^2)1' [L6) 

where |*i| = K^l^)] 1 / 2 . Solving the generalized eigenvalue problem with these Hamil- 
tonian and normalization matrix elements would yield improved upper bounds for the 
ground and low-lying excited states of the system. In the limit r — > oo the solutions 
would be exact. 
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In GFMC we can compute mixed expectation values such as 
(Mr)\ O \9 T j) 



(*i(r)|* r> i> 



(14) 



where the denominator involves just state i. Since the propagator commutes with the 
Hamiltonian, the desired matrix elements can be computed as: 



HiAr) 
NiAr) 



[H i:j H. 



l 



(15) 
(16) 



where we use expectation values computed from separately propagated ^(r) and tyj(r). 
For i = j these equations reduce to the standard GFMC calculation described above. 

As an example, Fig. El shows the computation of the energies of four | states in 7 Li. 
The lowest state has mainly [43] symmetry and can easily decay to the a+t channel; it has 
a large experimental width (918 keV) and its computed energy is slowly decreasing to the 
energy of the separated clusters. The remaining states are mainly of [421] symmetry and so 
are principally connected to the 6 Li+n channel. The second | state is experimentally just 
above the 6 Li+n threshold and has a small width (80 keV); its computed energy becomes 
constant with increasing r. The last two | states are not experimentally known, but the 
very slow decrease with r of the energy of the third state suggests that this state may also 
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Figure 3. GFMC energies of four | states in 7 Li versus imaginary time, r. 
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C IL2 result is preliminary. 



Figure 4. Energies of nuclear states computed with just the AV18 NN potential, and with 
the addition of the IL2 NNN potential, compared to experiment. 

be narrow. The off-diagonal overlaps Ni(t) are small and do not show signs of steadily 
increasing with increasing r. The solutions of generalized eigenvalue problems using the 
Eij(r) and iVy(r) are not significantly different from the Ea(r) shown in the figure. These 
results show that the (constrained) GFMC propagation largely retains the orthogonality 
of the starting Contrary to what might have been expected, the propagation of the 
higher states does not quickly collapse to the lowest state. 

4. ENERGIES OF NUCLEAR STATES 

Figure 0] compares energies computed with the AV18 (no V^) and AV18+IL2 Hamil- 
tonians to experimental values. The AV18+IL2 result shown for 12 C was made using 
a simplified \Pr and an approximate treatment of in the GFMC propagation; for 
these reasons it is marked preliminary. We see that using just a NN potential underbinds 
4 He by 4 MeV; this underbinding increases to 18 MeV for 12 C. The parameters of the 
Illinois-2 NNN potential were adjusted to reproduce the energies of 17 narrow states for 
3 < A < 8 [E]. As can be seen the potential provides an excellent overall reproduction of 
the energies of many states up to the ground state of 12 C; the RMS error in reproducing 
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Figure 5. Excitation energies of nuclear states. 



the experimental energies of the 39 states with width < 0.2 MeV is 0.7 MeV. 

Figure El shows excitation spectra of nuclei. One again sees a generally good repro- 
duction using AV18+IL2 of the experimental values. The NNN potential increases the 
predicted splittings of several spin-orbit partners, such as the 3 + ,2 + ,l + triplet in 6 Li and 
the first | and | levels in 7 Li which significantly improves the agreement with experi- 
ment. A variational study of 15 N also showed that a large part of the spin-orbit splitting 
in that nucleus is a result of the NNN potential [Ej. However the most dramatic effect 
of the NNN potential on levels appears in 10 B. The AV18 Hamiltonian results in two 1 + 
levels being below the 3 + state which is the experimental ground state; the AV18+IL2 
Hamiltonian gives the correct order. A similar inversion of levels occurs for the | and 
| states in 9 Be. These inversions are also due to increased spin-orbit strength from the 
NNN potential; in 1956 Kurath showed that the position of the 10 B(3 + ) state depends 
sensitively on spin-orbit strength [I15j. Another level ordering that is correct only with 
the NNN potential is shown in 10 Be. Here, as can be inferred from the experimental 
-B£7(2)'s to the ground state, the first 2 + level should have a large negative quadrupole 
moment and the second should have a positive quadrupole moment. This is the case for 
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the full AV18+IL2 Hamiltonian but without IL2 the energies of these states are reversed. 
5. OTHER RESULTS 

5.1. Comparison of GFMC and No-core Shell Model Calculations 

The no-core shell model (NCSM) is an alternative many-body method that uses real- 
istic NN and NNN potentials for light nuclei. It uses an expansion in a large harmonic- 
oscillator basis for all A nucleons. Effective two- and three-body interactions are con- 
structed from the bare u^-; bare Vyk can also be included in the three-body effective inter- 
actions. The construction of the effective interactions also gives a unique prescription for 
constructing effective operators from bare operators. The size of the harmonic-oscillator 
basis for a specified order of expansion grows rapidly with A. At present 16hQ can be 
used for A = 6 with just two-body effective interactions; this results in well- converged 
calculations. However only 8hQ can be used for A = 12, 16 and two-body effective in- 
teractions and these calculations are not fully converged. Calculations with three-body 
effective interactions are limited to 6hQ and are generally not converged. For a complete 
description of this method see [E] and references therein. 



-20 



-30 



-40 



-50 



& -60 




7 Li 

Argonne v 8 ' 
without & with TM'(99) 
NCSM & GFMC 



70 














: \ 




80 




\ AV8' 


AV8' 


AV8'+TM' 


AV8'+TM' 






1 NCSM 


GFMC 


NCSM 


GFMC 




c 











10 



B 



12 C 



Figure 6. Comparison of NCSM and GFMC energies for the AV8' and AV8'+TM' Hamil- 
tonians. 
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The NCSM is much faster than GFMC but because of these limits on the basis size 
appears to be less accurate than GFMC, especially when a NNN potential is being used. 
This is shown in Fig. El in which NCSM and GFMC calculations of several nuclei are 
compared. The left-hand two bars for each state show, respectively, NCSM and GFMC 
calculations for the AV8' potential [ 11] with no NNN potential; there is good agreement 
between both methods. The right-hand two bars shows results with the TM'(99) NNN 
potential [ E>] added to AV8'. In this case the NCSM results are significantly above 
those from GFMC; that is the NCSM appears to not get the full additional attraction 
provided by TM'(99). This is probably because of the limitation to only Qhfl when using 
the three-body effective interaction. 

5.2. Can Modern Nuclear Hamiltonians Tolerate a Bound Tetraneutron? 

An experimental claim of the existence of a bound tetraneutron cluster ( 4 n) was made 
recently [HZ|- This prompted our study [|Ej of the theoretical possibility of such a state; 
the conclusion was that a bound tetraneutron is strongly excluded by modern nuclear 
Hamiltonians. As a first step, negative-energy 4 n solutions using the AV18+IL2 model 
were searched for; GFMC calculations, using propagation to very large imaginary time 
(r = 1.6 MeV -1 ), produced only positive energies that steadily decreased as the RMS 
radius of the system increased. Adding artificial external wells to the Hamiltonian can, 
of course, produce negative energies. By varying the well depth, one can extrapolate 
to zero-well depth. Such calculations suggest that there might be a 4 n resonance near 
+2 MeV, but since the GFMC calculation with no external well shows no indication of 
stabilizing at that energy, the resonance, if it exists at all, must be very broad. 
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Figure 7. Predictions of Hamiltonians that have been modified to bind a tetraneutron. 
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To study the possibility that a minor modification of the AV18+IL2 model could result 
in binding of 4 n, a number of modifications to the AV18+IL2 model were made. In each 
case the modification was adjusted to bind 4 n with an energy of approximately —0.5 MeV; 
the consequences of this modification for other nuclei were then computed. Figure [7| shows 
some of these results. In the first such modification, the AVI' potential [ 18 j was used in 
just the partial wave (and AV18 in the other partial waves); this results in a 4 n energy 
of —0.52 MeV. However, as is shown by the bars labeled ^So in the figure, this results in 
substantial overbinding of 3 H, 4 He and all other nuclei that were computed. In addition 
it results in a bound dineutron; in fact the 4 n is not really bound as it can decay into two 
dineutrons. 

Modifications to the WWW or NNNN potentials, which are experimentally much less 
constrained than the AW potential, could be used to bind 4 n. A WWW potential that acts 
only in T = | triples would have the same effect on 4 n as one with no isospin dependence, 
but would have no effect on 3 H and 4 He because they contain only T — | triples. Such a 
potential was added to the AV18+IL2 Hamiltonian, and the coupling constant chosen to 
produce 4 n with ~ —0.5 MeV energy. It turns out that the coupling must be quite large 
to produce the minimally bound 4 n. 

This can be understood as follows. If a modified AW potential is used to bind 4 n, 
the pairs can sequentially come close enough to feel the attraction; this allows the four 
neutrons to be in a diffuse, large radius, distribution. However a NNN potential requires 
three neutrons to simultaneously be relatively close and thus the density of the system 
must be much higher. Indeed the RMS radii of the 4 n for the modified NNN potential is 
only 1.9 fm while the RMS radius for A n produced by the modified AW potential is 10 fm. 
The left panel of Fig. |H] shows the corresponding densities. The right panel compares the 




Figure 8. One- and two-body densities of 4 n bound with either a modified (circles) or 
T = 3/2 WWW potential (squares). The diamonds show twice the two-body density of 2 n. 
Two-body densities are shown projected into total spin or 1 channels. 
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pair distributions for the two different 4 n systems with that of the dineutron made by the 
modified NN potential. The 5 = pairs in the 4 n made by the modified NN potential are 
basically undisturbed 2 n pairs, while the modified NNN potential results in very different 
2 n pairs. 

The small radius of the 4 n bound by the modified NNN potential results in a kinetic 
energy that is an order of magnitude more than for the 4 n bound by the modified : 5o 
potential. This large kinetic energy must be overcome by a large NNN potential energy, 
and hence a large coupling constant is required. 

The very large coupling constant for the extra NNN term means that it has a large, 
even catastrophic, effect on any nuclear system in which it can act. This is shown in 
Fig- HI by the bars labeled V^; for example the binding energy of 6 Li is doubled and that 
of 6 He is tripled. As stated above, this potential has no effect on 3 H or 4 He. However 
the most dramatic result of this potential is that every investigated pure neutron system 
with A > 4 is extremely bound and in fact is the most stable "nucleus" of that A; the 6 n 
energy is off the scale of the figure at -650 MeV! 

A four-nucleon potential that acts in only T = 2 quadruples was also constructed to 
weakly bind 4 ; it has even more devastating consequences for nuclei in which it can act 
than the NNN potential. Thus an experimental observation of a bound tetraneutron 
would not be a minor perturbation to our understanding of nuclei; the experimental 
evidence for such a claim should be very strong before it is taken seriously. 



5.3. RMS radii of helium isotopes 

Recently a group at Argonne has measured the RMS charge radius of the radioactive 
nucleus 6 He (/3-decay halflife 0.8 sec.) with the remarkable accuracy of 0.7% [EH]. This 
prompted us to revisit our already published [ Ej proton RMS radius for 6 He. That 
calculation had been made with our standard propagation to r = 0.2 MeV -1 and had 
produced a value in excellent agreement with the new experimental value. We made 
new calculations to much larger r and found that there are very slow fluctuations of the 
RMS radius with r. These are shown in Fig. |U] which contains a GFMC propagation 
to t = 7.25MeV _1 for 6 He using the AV18+1L2 Hamiltonian. The upper panel shows 
the energy on a highly magnified scale; the results show just Monte Carlo statistical 
fluctuations with no long-term correlations. The standard error of the mean that is 
extracted from these numbers is thus reliable. However the RMS point proton radius 
shown in the lower panel has extremely long-range correlations making the error of the 
mean a useless number. We do not understand the origin of these fluctuations and have 
been unable to reduce them by using different importance sampling methods in the GFMC 
propagation. At the moment all we can do is report the average shown in the figure and 
use the extrema of the fluctuations as an estimate of the error; after folding in the proton 
and neutron RMS radii, this is reported in Table ^ as a charge radius. As can be seen 
our computed value is significantly above the experimental value. We have also made 
a large-r calculations of the RMS radii of 4 He and 8 He. The 4 He calculation shows no 
long-term fluctuations and that of 8 He shows substantially smaller fluctuations than for 
6 He. These results are also reported in the table. 
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Figure 9. GFMC propagation of 6 He energy (top panel) and RMS point proton radius 
(bottom panel). 

Table 1 

Computed and experimental RMS charge radii (in fm) of helium isotopes 





AV18+IL2 


Expt. 


4 He 


1.660(10) 


1.673(1) 


6 He 


2.15(7) 


2.054(14) 


8 He 


1.98(4) 





5.4. Center of Mass Effects on Spectroscopic Factors 

Spectroscopic factors are a measure of probability of finding subcluster structure in a 
nucleus. For the case of a specific A-l nuclear state in an A-body nucleus they are defined 
by computing the quasi-hole wave function for nucleon removal: 

X (r) = (^ A ^(J')xN(£ j )] J \a(r)\^A(J)) . (17) 
The spectroscopic factor is then 

S = J drr 2 \ X (r)\ 2 . (18) 

It is straightforward to compute S from shell-model wave functions. Conventional shell 
model calculations are done in a fixed center (FC) with harmonic oscillator wave functions. 
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Table 2 



Spectroscopic factors for Li(| ) — > He + p: + +p3, 2 + + pi, and 2 + +p3 computed 

2 2 2 

in a fixed center (FC) and with translational invariance (TI). 





7 T i 


l 

'r 2 \ 2 


<Sfc 


c 
o 


TI 




S P 




71 


0+ 


£ 


0+ 






1.95 1.95 


1.76 


1.84 


.56 


1. 


.65 


1.17 


7 

6 


1.95 2.5 


1.99 


2.16 


.56 


1. 


.63 


1.13 


1.13 


1.95 3.0 


2.22 


2.47 


.56 


1. 


.59 


1.04 


1.04 


Correlated 


2.27 


2.48 


.31 


0.56 


.33 


0.60 


1.07 


Correlated 


2.41 


2.50 


.32 


0.61 


.30 


0.54 


0.89 


VMC 


2.33 


2.48 






.36 


0.61 




Expt 


2.26 


"2.51" 






.42 


0.58 





Dieperink and de Forest showed how to convert such S to the desired translationally 
invariant (TI) ones. For the case of removing a single nucleon from the p shell, their 
formula becomes 

This increase of 5 arises from the center of mass acquiring a L = 1 component from the 
p-shell nucleons. This means that the nominally s-shell core nucleons have some p-wave 
component and thus contribute to the p-wave spectroscopic factors; s-wave spectroscopic 
factors are reduced to conserve the total number of nucleons. 

This theorem assumes that the same oscillator parameter is being used for s- and p- 
shell nucleons. However, this is not a realistic assumption; if the oscillator parameter is 
chosen to give a good RMS radius for the s-shell core ( 4 He) then very small RMS radii will 
be obtained for the p-shell nuclei. An example is shown in Table |2] which shows several 
calculations of spectroscopic factors for removing a proton from the ground state of 7 Li. 
The rows show various cases from an uncorrelated shell model to experimental values. 
The columns show the point RMS radii of 7 Li and the spectroscopic factors computed 
in a fixed center and with translationally-invariant wave functions. For each case the 
spectroscopic factors to the ground state (0 + ) of 6 He and the sum of the factors for the 
+ and 2 + states are shown. The final column shows the ratio of the TI to FC sums. 

The first three rows are for one-body harmonic-oscillator wave functions with no pair 
or triplet correlations; thus the FC results correspond to standard shell-model wave func- 
tions. The TI wave functions were made by expressing the oscillator wave functions in 
terms of — Rcm where is a nucleon coordinate and -Rcm is the center of mass of all 
A or A-l nucleons. The oscillator parameter for the s-shell gives a good RMS radius for 
4 He. This same parameter is used for the p-shell in the first line and results in very small 
7 Li RMS radii. In this case the full A/ (A — 1) ratio of TI to FC spectroscopic factors is 
observed, as is required by Eq. (|19|) . In the next two lines, the p-shell oscillator parameter 
is increased until reasonable 7 Li radii are achieved; this results in significantly reduced 



16 



TI/FC ratios. The two lines labeled "Correlated" report results from modified VMC 
wave functions which can be expressed in either the FC or TI. These have the structure 
of Eqn. (0|SJ) but also include one-body wave functions for the a-core nucleons. Here we 
see TI/FC ratios near 1 or even less than 1. The line labeled "VMC" shows the TI results 
for our full variational wave function; this wave function has no useful FC form. Based on 
these results it seems better to make no conversion of a FC spectroscopic factor computed 
with just a p-shell shell model, rather than increase it by A /{A — 1). We are continuing 
to study this. 

5.5. Mimicking Oxygen Isotopes as Neutron Drops 

Although GFMC is limited to light nuclei, it can be used to indirectly study large, 
neutron-rich, nuclei by computing neutron drops. These are collections of neutrons inter- 
acting via a realistic Hamiltonian, such as AV18+IL2, with the addition of an artificial 
external one-body well. The well provides the additional attraction necessary to produce 
a bound state of the N neutrons; it can be thought of as the average effect of the protons 
which are not explicitly included in the calculation. Note that the N neutrons occupy all 
shells from Os up; the well does not represent any omitted neutrons. 

Figure El shows results for one such study in which we attempt to mimic the neutron- 
rich isotopes of oxygen. Here we consider that the iV-neutron drop represents 8+N 0. The 
external well is of the Woods-Saxon form with parameters R = 3.0 fm, a = 1.1 fm, and 
V = —35.5 MeV; these were chosen in an attempt to get the correct separation energies 
for 17,18 0. Since the energy of the 8 protons is not included in the calculation, the energies 
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Figure 10. Energies of neutron drops compared to the corresponding oxygen isotope 
experimental energies. 
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of all the neutron drops are shifted (by 19 MeV) to match the 8-neutron drop with the 
experimental 16 O energy. The Woods-Saxon well binds Os and Op nucleons but does not 
have Od or Is bound states; thus the binding of neutrons in the sd shell arises from the 
combination of the well and the intrinsic AV18+IL2 potentials. 

We see that the energies of the neutron-rich oxygen isotopes up to 22 O are reasonably 
well reproduced by this model. However the 9-neutron drop has a | + ground state rather 
than the desired | + . It appears necessary to add a spin-orbit term to the well to achieve 
the correct ground state. Once a satisfactory model has been produced, the densities and 
other properties of these model oxygen isotopes can be studied. 

6. Conclusions 

Calculations with errors of only 1 to 2% of energies for nuclei from A = 6 to 12 for a 
given Hamiltonian are now possible. The AV18+IL2 Hamiltonian gives average binding- 
energy errors < 0.7 MeV for A = 3 — 12. A three-nucleon potential is required to obtain 
sufficient binding in the p-shell; it is also required to reproduce many experimental spin- 
orbit splittings and several level orderings. The GFMC and VMC wave functions can be 
used to study many other nuclear properties. 
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